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FOREWORD 


As  part  of  the  on-going  program  in  "Decision  Control  Models  in 
Operations  Research,"  Mr.  Charles  Mosier  has  refined  the  Power  Ap¬ 
proximation  for  computing  approximately  optimal  (s,S)  inventory 
policies.  Using  the  approach  adopted  by  Richard  Ehrhardt  (Technical 
Report  #7),  Mr.  Mosier  uses  regression  analysis  to  improve  the  accu¬ 
racy  of  an  analytically  derived  approximation.  Mr.  Mosier  refines 
the  Power  Approximation  by  constraining  the  regressions  to  provide  a 
policy  that  (1)  is  homogeneous  in  the  units  chosen  for  demand,  and 
(2)  has  reasonable  limiting  behavior  when  the  variance  of  demand  ap¬ 
proaches  zero.  The  improvements  are  obtained  with  only  modest  sacri¬ 
fices  in  total  cost  performance.  Other  related  reports  dealing  with 
this  research  program  are  given  on  the  following  pages. 


Harvey  M.  Wagner 
Principal  Investigator 


Richard  Ehrhardt 
Co-Principal  Investigator 
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algorithm  for  computing  (s,S)  inventory  policies.  The  approxima¬ 
tion  is  for  a  single-item,  periodic  review  model  with  set-up  cost, 
linear  holding  and  shortage  costs,  fixed  replenishment  lead  time, 
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to  derive  the  Power  Approximation  -  an  approximately  optimal  (s,S) 
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1.  INTRODUCTION  AND  SUMMARY 


1 . 1  The  Model 

We  consider  a  periodic  review,  single-item  inventory  system 
where  unfilled  demand  is  backlogged,  there  is  a  fixed  lead  time  L 
between  placement  and  delivery  of  an  order,  and  demands  during  re¬ 
view  periods  are  independently  and  identically  distributed,  having 

2 

mean  p  and  variance  a  .  Replenishment  costs  are  comprised  of  a 
setup  cost  K  and  a  unit  cost  c.  At  the  end  of  each  review  period 
a  cost  h  or  p  is  incurred  for  each  unit  on  hand  or  backlogged, 
respectively.  The  criterion  of  optimality  is  minimization  of  the 
undiscounted  expected  cost  per  period  over  an  infinite  horizon. 

Under  these  assumptions,  an  (s,S)  policy  is  optimal  (Iglehart 
[2]),  whenever  inventory  on  hand  or  backlogged  plus  on  order  y 
is  less  than  or  equal  to  s,  an  order  S-y  is  placed.  Algorith¬ 
mic  methods  for  computing  optimal  policies  are  available  (Veinott 
and  Wagner  [7]),  but  unfortunately,  the  computational  effort  re¬ 
quired  is  prohibitive  for  practical  implementation.  Furthermore, 
the  computation  of  an  optimal  policy  requires  the  complete  speci¬ 
fication  of  the  demand  distribution,  which  is  an  unrealistic  as¬ 
sumption  in  practical  settings.  Most  managers,  at  best,  have  im¬ 
precise  knowledge  of  the  mean  and  variance  of  the  demand  distribu¬ 
tion.  In  response  to  these  practical  limitations  in  available 
demand  information,  the  Power  Approximation  was  developed  by 
Ehrhardt  [1].  The  next  section  describes  this  approximation. 


i 

i 

i 


1 . 2  The  Power  Approximation 

The  Power  Approximation  computes  approximately  optimal  values 

for  (s,S)  using  only  the  mean  y  and  variance  a  of  demand. 

2  2 

Let  yL  =  (L+l)y  and  c?L  =  (L+l)a  .  The  algorithm  computes 

Dp  =  1.463y,364(K/h)>498a-138  ,  (1) 

2  =  {Dp/ [ ( 1  +  P/h)aL]}'5  ,  (2) 

and 

Sp  =  yL  +  ^l832(cj2/vi) '  18?( -220/z  +  1.142  -  2.866z)  .  (3) 

If  Dp/y  >  1.5,  let  S  =  Sp  +  Dp  and  s  =  Sp.  Otherwise,  compute 

s0  =  yL  +  vaL  ,  (4) 

where  v  is  the  solution  to 


-x^/2)/ /2irdx  =  p/(p+h)  . 


The  policy  parameters  are  given  by 


(5) 


s  =  minimum{Sp,So}  , 

S  =  minimum{s  +D  ,S  }  . 

r  r  y 

If  demands  are  integer  valued,  s  ,  D  ,  and  S  are  rounded  to 

r  r  u 

the  nearest  integer. 


1.3  Motivation  for  This  Study 

Ehrhardt  [1]  discussed  how  theoretical  considerations  lead  to 


the  following  form  for  an  approximately  optimal  policy. 


D  5  (2Ky/h) 


(6) 


.5 


and 


s  =  uL  +  aLG(z)  ,  (7) 

where 

z  £  D/[(l  +  p/h)oL]  ,  (8) 

and  G(*)  Is  a  function  that  depends  upon  the  demand  distribution. 
The  Power  Approximation  was  derived  by  using  these  expressions  as 
the  bases  of  regression  analyses,  where  known  optimal  policies  are 
used  as  data.  The  regressions  yielded  parameter  values  in  func¬ 
tions  having  the  general  forms  (6),  (7),  and  (8).  The  resulting 
expressions  (1),  (2),  and  (3)  are  further  modified  by  (4)  and  (5) 
according  to  the  empirical  modification  of  Wagner  [7],  In  this 
study,  we  focus  only  on  modifying  (1),  (2),  and  (3)  to  overcome 
several  deficiencies  in  the  original  analysis. 

The  Power  Approximation  has  three  deficiencies.  First,  Dp 
in  (1)  converges  to  zero  as  the  variance  of  demand  goes  to  zero, 
whereas  we  would  prefer  that  Dp  converge  to  the  Wilson  Lot  Size 

s  (2Ky/h)'5.  Second,  the  approximation  is  not  homogeneous  with 
w 

respect  to  the  scaling  of  demand.  Third,  the  numerical  analysis 
conducted  for  fitting  the  approximation  erred  In  a  minor  detail. 
These  deficiencies  are  discussed  in  the  following  three  sections. 
1.3.1  Limit  of  D 

The  condition  that  11m  Dp  *  Dw  ■  (2Ky/h)‘®  as  o2  0, 
does  not  hold  in  (1)  for  the  Power  Approximation.  In  fact, 
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2 

lim  Dp  =  0  as  a  ■*  0.  This  deficiency  is  remedied  by  using  the 
following  form 

Dp  =  au(1’B)(K/h)eO  +  a[/y2)Y  .  (9) 

Our  analysis,  as  described  in  the  following  sections,  yields  the 
result, 

Dp  =  1 . 30y ‘ 494 ( K/h ) ‘ 506 ( 1  +  ajVy2)*116  .  (10) 

Hence , 

lim  D  =  1.30y'494(K/h)-506  ,  (11) 

o+O 

which  is  very  close  to  the  Wilson  Lot  Size 

Dw  *  1.414y-5(K/h)-5  .  (12) 

1.3.2  Non- homogeneity  Deficiency 

If  demand  is  rescaled  by  a  factor,  say  k,  then  s  and  D 
should  be  transformed  similarly.  That  is,  if  y  =  ky  and 
a  *  ka,  then  we  should  have  that  D  =  kD  and  s  =  ks.  The  Power 
Approximation,  (1),  (2),  and  (3),  however,  are  not  homogeneous. 

We  remedy  this  situation  in  the  present  analysis  as  follows. 

In  [1],  the  regression  for  Dp  used  the  model 

Dp  =  ay°(K/h)6(oL)Y  .  (13) 

In  the  current  analysis,  expression  (9)  is  the  form  used;  thus,  we 
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have  replaced  (cr^)Y  with  (1  +  o2/y2)Y.  The  1  added  to  the 

A  2 

last  term  ensures  that  Dp  is  not  zero  in  the  limit  and  the  y 
divisor  forces  the  last  term  to  remain  constant  on  a  rescaling  of 
demand.  Requiring  that  a  +  8  =  1  forces  homogeneity,  since  if 
we  rescale  demand  by  k,  then  the  new  mean  is  k  times  the  old 
mean  and  the  new  holding  cost  is  equal  to  the  old  holding  cost 
divided  by  k. 

Homogeneity  was  preserved  in  the  fit  for  s  by  requiring  a 
fit  of  the  form  in  (7),  where  6(*)  is  dimensionless.  This  dif¬ 
fers  from  the  original  Power  Approximation,  where  was 
replaced  with  a  different  function  of  demand  parameters. 

1.3.3  Error  in  Numerical  Analysis 

In  the  original  numerical  analysis  [1],  three  candidates  for 
z  in  expression  (3)  were  considered.  They  were 

»  <Dp/((p/h)oL))-5  ,  (14) 

z2  =  ^Dp  +  +  °2Ai2))/(0  +  P/h)oL))‘5  ,  (15) 

and 

z3  s  (yW  +  P/h)\))‘5  •  06) 

Expression  (16)  was  chosen  over  (14)  and  (15)  because  It  provided 
the  best  numerical  fit  to  optimal  policy  data.  Inadvertantly, 
expression  (15)  was  computed  as 

z2  -  ((Dp  +  .5(1  +  o2/y2))/(l  +  p/h)oL)'5  (17) 

In  the  original  analysis.  In  the  current  analysis,  the  correct 
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expression  is  considered.  It  does  not,  however,  provide  the  best 
numerical  fit  with  the  data. 

2.  EXPERIMENTAL  DESIGN 

The  data  for  this  study  are  the  same  as  those  in  the  original 
analysis  of  the  Power  Approximation  [1].  An  experimental  grid  of 
288  parameter  settings  is  specified  and  presented  in  Table  1. 

Three  types  of  demand  distributions  are  examined,  namely, 
Poisson  and  Negative  Binomial  with  variance-to-mean  ratios  of  3 
and  9.  Each  demand  distribution  is  evaluated  with  four  mean 
values,  2,  4,  8,  and  16.  Lead  time  has  three  values,  0,  2, 
and  4.  Since  the  cost  function  is  linear  in  the  parameters  K, 
p,  and  h,  the  unit  holding  cost  value  can  be  normalized  at 
unity.  The  unit  penalty  costs  are  4,  9,  24,  and  99,  and  the 
set  up  cost  values  are  32  and  64.  The  unit  replenishment  cost 
c  need  not  be  specified,  since  it  does  not  affect  the  computation 
of  an  optimal  policy  for  an  undiscounted,  infinite  horizon  model 
with  complete  backlogging.  All  combinations  of  the  parameter 
settings  are  included  in  the  grid,  yielding  288  points. 

The  optimal  policy  for  each  of  the  288  settings  is  calculated 
by  the  algorithm  of  Veinott  and  Wagner  [7],  implemented  in  the 
software  written  by  Kaufman  [4].  The  resulting  288  values  of  s 
and  S  are  the  data  utilized  in  the  least  squares  regressions  and 
the  evaluation  of  results.  The  regression  analysis  is  performed 
using  the  Statistical  Package  for  the  Social  Slcences,  second 
edition  (1975). 
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TABLE  1 

System  Parameters 


Factor 

;;  I 

Levels  i 

Demand  Distribution 

*  Negative  Binomial 

{Variance-to-Mean  Ratio  =  9) 

Negative  Binomial 
(Variance-to-Mean  Ratio  =  3) 

Poisson 

(Variance-to-Mean  Ratio  =  1) 

Mean  Demand  (y) 

2,  4,  8,  16 

Replenishment  Lead 
Time  (L) 

0,  2,  4 

Replenishment  Setup 
Cost  (K) 

32,  64 

Unit  Penalty  Cost 
(p) 

4,  9,  24,  99 

Unit  Holding  Cost 

00 

1 

No.  of 
Levels 


3.  APPROXIMATIONS  FOR  D 


Several  mathematical  forms  were  tested  in  seeking  a  new 
approximation  for  D.  We  considered  a  linear  model 

D  =  (2Ky/h) "5  +  [A  +  Bo/y  +  Co2/y2]  ,  (18) 


and  a  multiplicative  model 


D  =  aya(K/h)e(l  +  a2/y2)Y(p/h)6  .  (19) 

The  best  fits  were  obtained  with  the  multiplicative  form.  Initial 
fits  showed  the  variable  p/h  to  be  insignificant,  so  we  arrived 
at  the  final  model 

D  =  aya(K/h)e(l  +  a2/y2)Y  ,  (20) 

where  we  require  a  +  B  =  1  (see  Sec.  1.3.2). 

Taking  logarithms  in  (20),  we  have  the  linear  expression 

£nDp  =  In  a  +  a  in\i  +  &ln  (K/h)  +  y£n  (1  +  a2/y2)  .  (21) 


We  let  a  =  1  -  3,  yielding 


£nDp  =  a  +  (1  -  8 )  y  +  8  &i  ( K/h )  +  y  ( 1  +  cr£ A*  ) 
The  final  form  for  fitting  is 


In  D  -  £ny  -  In  a  +  8  [In  (K/h)  -  £ny]  +  yin  (1  +  o^/y  / 


The  result  is 


0  *  1.30y,494(K/h),506(l  +  a2/y2)’116 


with  R  =  .982. 
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4.  APPROXIMATIONS  FOR  s„ 

P 

The  general  form  for  s  (from  [1])  is 

sp  =  ayL  +  aL(A/z^  +  B  +  Cz^ )  . 
Three  alternative  forms  of  z  were  considered: 

=  C0p/({p/h)aL)]-5  , 

[Dp  +  .5(y  +  cr2/y2)" |'5 

Z2  =  [  (1  +  p/h)aL  J 

z3  =  [Dp /( (1  +  p/h)aL)]'5  . 
Each  satisfies  the  assumptions  of  Roberts  [6]. 
Model  1  (zj). 


z  =  [Dp/((p/h)oL)]'5  . 

Using  least  squares  regression,  we  have 

Sp  =  .973yL  +  aL(.183/z  +  1.063  -  2.192z)  , 
with  R2  =  .997. 


Model  2  ( z 2 )  . 


Let 


Z  *  [(Dp  +  • 5 (y  +  o2/y2))/((l  +  p/h)oL)]*5  . 


(25) 

(26) 

(27) 

(28) 


(29) 


(30) 


(31) 


>  itniTiWaAiM  f  i  ■ 


Then  least  squares  regression  yields 


sp  =  .929yL  +  ol(.190/z  +  1.263  -  2.410z)  ,  (32) 

with  R2  *  .998. 


Model  3  Ug)- 


Let 

z  =  [Dp/( (1  +  p/h)oL)]*5  .  (33) 

Least  square  produces 

Sp  =  .977yL  +  ol(.171/z  +  1.174  -  2.652z)  +  .36  ,  (34) 

which  we  call  Model  3A;  this  model  has  R2  =  .997.  We  also  consider 
a  fourth  rounded-down  version,  denoted  Model  3B, 

sp  =  .977yL  +  0^.171/z  +  1.174  -  2.652z)  .  (35) 

4.1  Cost  Comparison  and  Choice  of  the  Final  Model 

The  288  settings  used  in  the  regression  analyses  are  examined 
first.  The  expected  total  cost  per  period  is  calculated  for  each 
item  when  controlled  with  each  of  the  five  policies:  the  optimal 
policy  and  the  four  policy  generating  models.  All  policy  generating 
models  evaluate  Dp  using  expression  (24),  and  sp  using  expres¬ 
sions  (30),  (32),  (34),  and  (35)  for  Models  1,  2,  3A,  and  3B, 
respectively. 

Let  C(m),  for  m  ■  1,  2,  3A,  3B,  and  C*  be  the  expected 
total  cost  per  period  for  an  item  when  controlled  using  a  policy 
approximation  (each  of  the  four  models)  and  the  optimal  policy. 
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respectively.  Our  criterion  of  performance  is 

Ap  =  { [C (m)  -  C*]/C}*100%  ,  (36) 

which  is  the  percentage  by  which  the  expected  cost  exceeds  the 
optimal  cost.  The  results  for  the  238  settings  are  summarized 
in  Table  2,  which  lists  the  number  of  cases  having  values  of  Ap 
in  various  ranges. 

Ap  is  displayed  for  the  Power  Approximation  fcr  comparitive 
purposes.  Table  3  gives  cumulative  percentages. 

Table  4  examines  the  characteristics  of  the  settings  which 
seem  to  be  "outliers,"  that  is,  settings  with  Ap  at  least  3%. 

The  best  of  the  four  models  seems  to  be  Model  1.  The  average 
error  of  Model  1  is  lowest.  The  number  of  items  of  low  accuracy, 
for  example,  an  error  greater  than  3.0%,  is  smallest.  Finally, 
for  a  majority  of  the  outliers,  the  error  for  Model  1  is  less  than 
or  equal  to  that  for  the  other  models. 

5.  REVISED  POWER  APPROXIMATION 


The  algorithm  for  an  approximately  optimal  (s,S)  policy  is  as 
follows.  Compute 

Dp  =  1.30M’494(K/h),506(l  +  a2/y2)-116  (37) 

and 


sp  =  .973yL  +  ol(.183/z  +  1.063  -  2.192z)  ,  (38) 


where 


z  -  [Dp/((p/h)oL]'5  . 


(39) 
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TABLE  2 

ERROR  FREQUENCIES 


Model 

AP 

1 

2 

3A 

3B 

[O.O. 0.1) 

151 

118 

117 

108 

116 

[0.1, 0.5) 

102 

108 

107 

110 

105 

[0.5, 1.0) 

21 

25 

31 

31 

32 

[1.0, 2.0) 

11 

18 

11 

16 

15 

[2. 0,3.0) 

3 

11 

13 

14 

12 

[3. 0,4.0) 

0 

4 

4 

6 

2 

[4. 0,5.0) 

0 

3 

3 

2 

4 

[5. 0,6.0) 

0 

1 

2 

1 

2 

Avg.  Error 

0.35% 

0.469% 

0.438% 

0.507% 

0.484% 

TABLE 

3 

CUMULATIVE  FREQUENCY 

(%)  OF  ERRORS 

Model 

.  ap  . 

Power 

Approx. 

1 

2 

3A 

3B 

[0. 0,0.1) 

52% 

41% 

40.6% 

37.5% 

40.3% 

[0.1,0. 5) 

87.8% 

78.5% 

77.8% 

75.7% 

76.7% 

[0.5, 1.0) 

95% 

87.2% 

88.5% 

86.5% 

86.5% 

[1.0, 2.0) 

99% 

93.4% 

92.4% 

92.0% 

93.1% 

[2. 0,3.0) 

100% 

97.2% 

96.9% 

96.9% 

97.2% 

[3. 0,4.0) 

100% 

98.6% 

98.3% 

99.0% 

97.9% 

[4. 0,5.0) 

100% 

99.7% 

99.3% 

99.7% 

99.3% 

[5. 0,6.0) 

100% 

100% 

100% 

100% 

100% 
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( 

r 

j  TABLE  4 

OUTLIERS 


*Not  an  outlier  for  this  particular  model. 


Otherwise 


If  D_/]j  >  1.5,  then  let  S  =  S_  +  and  s  =  s  . 
p  P  P  P 

compute 

S0  =  yL  +  V0L  • 
where  v  is  the  solution  to 


-x2/2)//2tF  dx  =  p/(p+h)  . 


The  policy  parameters  are 


and 


s  =  minimum  |Sp,Sg}  , 


(40) 


(41) 


S  =  minimum  {sp+Dp,S0}  . 

If  demands  are  integer  valued,  then  sp,Dp,  and  sq  are  rounded 
to  the  nearest  integer. 

5.1  Robustness  of  Revised  Power  Approximation 

Here  we  examine  parameter  settings  other  than  those  used  in 
deriving  the  approximation.  Table  5  describes  the  settings  for  a 
full -factorial  design  with  32  cases,  having  both  interpolated  and 
extrapolated  values  for  the  parameters.  With  the  exception  of  the 
unit  holding  cost,  each  parameter  Is  set  at  two  new  values:  one 
Is  an  interpolation  between  levels  in  Table  1,  and  the  other  Is  an 
extrapolation  beyond  previous  values  in  Table  1.  Table  6  lists 
the  frequencies  of  Ap  for  the  32  cases  and  also  lists  the  re¬ 
sults  of  the  original  Power  Approximation.  The  costs  for  both  ap¬ 
proximations  are  only  slightly  higher  than  optimal;  the  costs  for 
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TABLE  5 

INTERPOLATED  AND  EXTRAPOLATED  VALUES 


Factor 

Levels 

No.  of 
Levels 

Demand  Distribution 

Negative  Binomial 
Variance-to-Mean  Ratio  =  5 

2 

Negative  Binomial 
Variance-to-Mean  Ratio  =  15 

Mean  Demand  (p) 

0.5,  7.0 

2 

Replenishment  Lead 
Time  (L) 

1,  6 

2 

Replenishment  Setup 
Time  (K) 

16,  43 

2 

Unit  Penalty  Cost 
(P) 

49,  132 

2 

Unit  Holding  Cost 
(h) 

1 

1 

TABLE  6 

ERROR  FREQUENCY  AND  CUMULATIVE  PERCENTAGE 
FOR  INTERPOLATED  AND  EXTRAPOLATED  PARAMETER  SETTINGS 


AP 

Power  Approximation 

Model  1 

Frequency 

Percentage 

Frequency 

Percentage 

[0.0%, 0.2%) 

12 

38% 

8 

25% 

[0.2%, 0.6%) 

S 

61% 

4 

37.5% 

[0.6%, 1.0%) 

4 

75% 

2 

43.8% 

[1.0%, 2.0%) 

3 

84% 

5 

59.4% 

[2.0%, 3.0%) 

3 

94% 

2 

65.6% 

[3.0%, 4.0%) 

1 

97% 

6 

84.4% 

[5.0%, 11.0%) 

1 

100% 

5 

100% 

£  ^  ^ 
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the  new  approximation  are  slightly  higher  than  those  for  the  original 
Power  Approximation.  Of  the  eleven  settings  with  values  of  Ap 
greater  than  3.0%,  two  have  two  parameters  with  extrapolated 
values,  four  have  four  parameters  with  extrapolated  values,  and  one 
has  five  parameters  with  extrapolated  values. 

We  also  examine  the  extreme  extrapolations  for  individual 
parameter  values.  A  base  case  is  chosen  for  comparison.  The  param¬ 
eter  settings  of  the  base  case  are  near  the  midpoints  of  the  ranges 
used  in  the  288  settings  (negative  binomial  demand,  a  /y  =  5, 
y  =  9,  L=2,  h=l,  p  =  49,  and  K  =  48).  The  value  of  the 
variance-to-mean  ratio,  lead  time,  and  penalty  cost  were  extrapo¬ 
lated  to  20,  10,  and  1 99--sl i ghtly  more  than  double  the  largest 
values  used  in  the  288  settings.  Table  7  lists  the  parameter 
settings  and  the  resulting  Ap  for  each  of  the  extrapolations  for 
both  the  Power  Approximation  and  Model  1. 

We  see  in  all  cases  that  the  Power  Approximation  yields  total 
costs  within  1%  of  optimal,  and  Model  1  has  only  one  case  of 
slightly  higher  costs. 


6.  CONCLUSIONS 

We  have  derived  an  approximately  optimal  policy  that  is  easily 
computed,  requires  only  the  mean  and  variance  of  the  demand  distri¬ 
bution,  and  provides  a  good  approximation  to  optimality  over  a  wide 
range  of  parameter  settings.  The  approximation  is  accurate  when  com¬ 
pared  to  optimal  and  is  only  slightly  more  costly  than  the  original 
Power  Approximation  due  to  the  added  constraints  ensuring  homogeneity 
and  the  proper  limit  for  approximation  of  D. 
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TABLE  7 

SINGLE  PARAMETER  EXTRAPOLATIONS 

Base  Case:  Negative  Binomial  Demand 

(Variance-to-Mean  Ratio  =  5, 
y  =  9,  L  =  2,  p  =  49,  K  =  48) 


Power 

Extrapolated  Value  Approximation  Model  1 


Variance-to-Mean 

Ratio  =  20 

0.0% 

1.06% 

y  =  20 

0.10% 

0.01% 

30 

0.21% 

0.14% 

40 

0.18% 

0.16% 

K  =  20 

0.11% 

0.00% 

15 

0.28% 

0.16% 

9 

0.63% 

0.43% 

P  =  132 

0.15% 

0.02% 

199 

0.50% 

0.18% 

L  =  10 


0.02% 


0.03% 
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